f(T) gravity mimicking dynamical dark energy. Background and perturbation analysis 
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We investigate f(T) cosmology in both the background, as well as in the perturbation level, and we 
present the general formalism for reconstructing the equivalent one-parameter family of f(T) models 
for any given dynamical dark energy scenario. Despite the completely indistinguishable background 
behavior, the perturbations break this degeneracy and the growth histories of all these models differ 
from one another. As an application we reconstruct the f(T) equivalent for quintessence, and we 
show that the deviation of the matter overdensity evolution is strong for small scales and weak for 
large scales, while it is negligible for large redshifts. 
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I. INTRODUCTION 

Cosmological data from a wide range of sources in- 
cluding type la supernovae [IH3, the cosmic microwave 
background jjK baryon acoustic oscillations [a, M , cluster 
gas fractions 0, Q and gamma ray bursts [3jlo| seem to 
indicate that at least 70% of the energy density in the 
universe is in the form of an exotic, negative-pressure 
component, which drives the universe acceleration. Al- 
though the simplest way to explain this behavior is the 
consideration of a cosmological constant the known 
fine-tuning problem led to the dark energy paradigm. 
The dynamical nature of dark energy, at least in an ef- 
fective level, can originate from a variable cosmological 
"constant" [12j , or from various fields, such is a canonical 
scalar field (quintessence) (l3l - [l7l |. a phantom field, that 
is a scalar field with a negative sign of the kinetic term 
[l8j . the combination of quintessence and phantom in a 
unified model named quintom or from k-essence [2(| 
and unparticles [2l|. On the other hand, the dynamical 
dark energy can effectively be described by modifying 
gra vity itself, using a functions of the curvature scalar 
|22| , of the Gauss-Bonnet invariant [23| , or of the square 
of the Weyl tensor [24||, hi gher derivatives in the action 
(25| , braneworld extensions |2fj| , string- inspired construc- 
tions holographic properties (28|, UV modifications 
[H etc. 

An interesting alternative to General Relativity is the 
so-called f(T) gravity, which has recently received con- 
siderable attention as a possible explanation of the late- 
time acceleration |30h33 |. It is based on the old idea 
of the "teleparallel" equivalence of General Relativity 
(TEGR) |33l, |34|. which, instead of using the Reimann- 
Cartan space with curvature but no torsion defined via 



the Levi-Civita connection, uses the Weitzenbock con- 
nection that has no curvature but only torsion. The dy- 
namical objects in such a framework are the four linearly 
independent vierbeins. The advantage of this framework 
is that the torsion tensor is formed solely from products 
of first derivatives of the tetrad. As described in [34| , the 
Lagrangian density, T, can then be constructed from this 
torsion tensor under the assumptions of invariance under 
general coordinate transformations, global Lorentz trans- 
formations, and the parity operation, along with requir- 
ing the Lagrangian density to be second order in the tor- 
sion tensor. However, instead of using the tensor scalar 
T the authors of [33, |32| generalized the above formalism 
to a modified f(T) version, thus making the Lagrangian 
density a function of T, similar to the well-known exten- 
sion of f(R) Einstein- Hilbert action. In comparison with 
f(R) gravity, whose fourth-order equations may lead to 
pathologies, f(T) gravity has the significant advantage of 
producing field equations which are at most second or- 
der in field derivatives. This feature has led to a rapidly 
increasing interest in the literature, and apart from ob- 
taining acceleration [3lH32| one can reconstruct a va riety 
of cosmological evolutions |35l.[36j. add a scalar field [37l |. 
examine the conformal transformations [38[, use obser- 
vational data in order to constrain the model parame- 
ters |39| , examine the dynamical behavior of the scenario 
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4C| and the possibility of the phantom divide crossing 
41 , |42j , and proceed beyond the background evolution, 
investigating the vacuum and matter perturbations [43[ . 

In this paper we are interested in constructing f(T) 
scenarios, that exhibit the same background behavior 
with any given cosmological model. However, we ad- 
ditionally investigate the perturbation evolution, and in 
particular we examine the growth of matter overdensity, 
since it distinguishes between the dynamical dark energy 
scenario and its equivalent family of f(T) models. 

The layout of this paper is as follows. In Section [II] 
we briefly review the cosmology of f(T) gravity, both 
at the background and linearized regimes. In Section 
IIIII we present the general formalism of reconstructing 
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the equivalent f(T) models of any given dynamical dark 
energy scenario, as well as the formalism of using per- 
turbations in order to distinguish them. In Section II VI 
we apply the obtained formalism in the quintessence sce- 
nario. Finally, Section fVl summarizes our results. 

II. fiT) GRAVITY AND COSMOLOGY 

In this section we present f(T) gravity, we provide 
the background cosmological equations in a universe gov- 
erned by such a modified gravitational sector, and we 
give the first order perturbed equations. Throughout 
the work we consider a flat Friedmann- Robertson- Walker 
(FRW) background geometry with metric 

ds 2 = dt 2 - a 2 {t) 5 lJ dx l dx j , (1) 

where a(t) is the scale factor. In this manuscript our 
notation is as follows: Greek indices run over all 

coordinate space-time 0, 1, 2, 3, lower case Latin indices 
(from the middle of the alphabet) ... run over spatial 
coordinates 1, 2, 3, capital Latin indices A, B,... run over 
the tangent space-time 0, 1, 2, 3, and lower case Latin 
indices (from the beginning of the alphabet) a, will 
run over the tangent space spatial coordinates 1, 2, 3. 

A. f(T) gravity 

Let us present f(T) gravity. As stated in the Intro- 
duction, the dynamical variable of the old "teleparallel" 
gravity, as well as its f(T) extension, is the vierbein field 
eA(x^), where capital as well as Greek indices take the 
values 0, 1, 2, 3. This forms an orthonormal basis for the 
tangent space at each point x^ of the manifold, that 
is e A ■ e B = rj AB , where r] A B = diag(l,—l,—l,—l). 
Furthermore, the vector can be analyzed with the 
use of its components e A in a coordinate basis, that is 
e A = e^<9 M . 

In such a framework, the metric tensor is obtained from 
the dual vierbein as 

9nv (x) = ijab e-l (x) ef (x) . (2) 

Contrary to General Relativity, which uses the torsion- 
less Levi-Civita connection, in the present formalism 
ones uses the curvatureless Weitzenbock connection [44| , 
whose torsion tensor reads 

r^ = r^-r^ = el(9^-^). (3) 

Moreover, the contorsion tensor, which equals the differ- 
ence between Weitzenbock and Levi-Civita connections, 
is defined as 

K^ p = -l(T^ p -T v »-T p » v ). (4) 
Finally, it proves useful to define 

S p = \ (K% + 5% T au a 6; T°*) . (5) 



Note the antisymmetric relations T x fll , = — T A „ M and 
Sp^ = —Sp vtl , as can be easily verified. Using these 
quantities one can define the so called "teleparallel La- 
grangian" as 0, [H, Si 

L T ^V" T V (6) 

In summary, in the present formalism, all the informa- 
tion concerning the gravitational field is included in the 
torsion tensor , and the teleparallel Lag rangian Lj* 
gives rise to the dynamical equations for the vierbein, 
which imply the Einstein equations for the metric. 

From the above discussion one can deduce that the 
teleparallel Lagrangian arises from the torsion tensor, 
similar to the way the curvature scalar arises from the 
curvature (Riemann) tensor. Thus, one can simplify the 
notation by replacing the symbol Lt by the symbol T, 
which is the torsion scalar [32| • 

While in teleparallel gravity the action is constructed 
by the teleparallel Lagrangian Lt = T, the idea of /(T) 
gravity is to generalize T to a function T + f(T), which is 
similar in spirit to the generalization of the Ricci scalar R 
in the Einstein-Hilbert action to a function f(R). In par- 
ticular, the action in a universe governed by f(T) gravity 
reads: 

I=j^ J d 4 xe[T + f(T) + L m ], (7) 

where e = det(e^) = y/—g and L m stands for the matter 
Lagrangian. We mention here that since the Ricci scalar 
R and the torsion scalar T differ only by a total derivative 
(47| . in the case where f(T) is a constant (which will 
play the role of a cosmological constant) the action ([7]) 
is equivalent to General Relativity with a cosmological 
constant. 

Finally, we mention that throughout this work we use 
the common choice for the form of the vierbien, namely 

= diag(l, a, a, a). (8) 

It can be easily found that the family of vierbiens related 
to (JSJ) through global Lorentz transformations, lead to 
the same equations of motion. Note however that, as it 
was shown in [48[, f(T) gravity does not preserve local 
Lorentz invariance. Thus, one should in principle study 
the cosmological consequences of a more general vierbien 
ansatz. Definitely the subject needs further investigation 
and it is left for a future work. 



B. Background f(T) cosmology 

Let us now present the background cosmological equa- 
tions in a universe governed by f(T) gravity. Variation 
of the action ([7]) with respect to the vierbein gives the 
equations of motion 
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-%(eS A n[i + /'CO] - e A T% x S p ^[l + f(T)] + S A ^(T)f"(T) - -e A [T + f(T)] 



em 

4nGe p A T p \ 



(9) 



r 



where a prime denotes the derivative with respect to T 



and the mixed indices are used as in Sa^ 



Note that the tensor T p v on the right-hand side is the 
usual energy-momentum tensor, in which we have added 
an overset label in order to avoid confusion with the tor- 
sion tensor. 

If we assume the background to be a perfect fluid, then 
the energy momentum tensor takes the form 



em 



V9nv - (p+pjUftU,,, (10) 

where is the fluid four-velocity. Note that we are 
following the conventions of |47j, but with an opposite 
signature metric. Thus, one can see that the equations 
^ lead to the background (Friedmann) equations 



H 2 



H = 



8ttG 



/CO 
; r "' 6 
4nG(p m + 



-/>„ 



- 2f'{T)H 2 (11) 

Pm) 



(12) 

1 + f'(T) — 12H 2 f"(T) ' v ' 

In these expressions we have introduced the Hubble pa- 
rameter H = a/ a, where a dot denotes a derivative with 
respect to coordinate time t. Moreover, p m and p m stand 
respectively for the energy density and pressure of the 
matter content of the universe, with equation-of-state pa- 
rameter w m = Pm/pm- Finally, we have employed the 
very useful relation 

T = -6H 2 , (13) 

which straightforwardly arises from evaluation of ^ for 
the unperturbed metric. 

Observing the form of the first Friedmann equation 
pip , and comparing to the standard form, we deduce 
that the second and third terms on the right hand side 
constitute effectively the dark energy sector, which in 
general presents a dynamical behavior. In particular, one 
can define the dynamical dark energy (DDE) density as 
HI: 



Pdde 



8ttG 



/CO 
6 



2f(T)H 2 



(14) 



(15) 



while its equation-of-state parameter reads: 

f IT — f'(T) + 2Tf"(T) 

w = . 

[1 + f'(T) + 2T/"(T)] [f/T - 2/'(T)] 

Thus, in principle, any dynamical dark energy scenario, 
with a given pdde or a given w, has its f(T) equivalent, 
and the corresponding reconstruction procedure will be 
described in section Hill Lastly, note that General Rela- 
tivity is recovered by setting f(T) to a constant (which 
will play the role of a cosmological constant), as expected. 
C. Linear Perturbations in f(T) gravity 

We now recall the first order perturbations of f(T) 
gravity. We refer the reader to [43| for full details of the 



calculation, summarizing only the relevant results in this 
section. The perturbed metric in Newtonian gauge with 
signature (H ) is written as 



ds 2 



(1 + 2i>)dt 2 - a 2 (l - 24>)5 ij dx i dx j , 



(16) 



where the scalars <j> and ip are functions of x and t. Writ- 
ing the vierbein as 



e A = B A + t A 



(17) 



where e A is the unperturbed vierbein, for the above met- 
ric we obtain 



6?, e?, =8> i$ = 88 e£ 



SH 



(18) 



8X 



tl = W t a » = S*acb t% = -8^ ^ = ^0,(19) 

with indicial notation as stated at the beginning of sec- 
tion [ill Moreover, unless otherwise indicated, subscripts 
zero and one will generally denote respectively zeroth and 
linear order values of quantities. We have also incorpo- 
rated the additional simplifying assumption that the per- 
turbations t A are diagonal. With these perturbations the 
determinant becomes 

e = det(e^) = a 3 (l + ip -3<f>). (20) 

The perturbations of the energy-momentum tensor are 
expressed as [43[ 

em 

(5 T ° = S Pm (21) 

em 

ST a = a - 2 (p m +p m )(~d a Su) (22) 

em 

ST a ° - { Pm + PmKdaSu) (23) 



ST a = S ab 5p. 



(24) 



Pressureless matter implies that (j> = ip, a relation that 
simplifies the calculations significantly. Additionally, it 
proves convenient to transform to Fourier space, where 
any first order quantity y is expanded as 



d 3 k 
(2tt)1 



(25) 



In what follows, the k-subscripts in a quantity denotes 
its Fourier transformation. 

Finally, from the above one can result to the following 
evolution equation for the metric perturbation (43| : 



^ + r(/)0 fc + w 2 (/)^ fe = o, 

where uj 2 and T are respectively given by 



(26) 
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3H S 



H 



Hf - ( 36ff 4 - A8H 2 H) f" + UAH 4 Hf 



1 + /' - 12H 2 f" 



(27) 



r(/) 



AH 



l + f- [12H 2 + 9Hj f" + 3QH 2 Hf 
l + f - 12H 2 f> 



(28) 



with primes denoting derivatives with respect to T. 



D. Growth of perturbations 

In order to study the growth of perturbations in a 
matter-only universe (p m = 6p m = 0) we define as usual 
the overdensity 5 as 



(29) 



Thus, the relativistic version of the Poisson equation in 
f(T) gravity reads [1| 

3H (1 + f - 12H 2 f) <j> k 
+ [(3H 2 + k 2 /a 2 ) (1 + /') - 36ff 4 /"] <j> k 
+ AnGp5 k = Q. (30) 

Equations (f2l?)) and ([3U]) can be used to evolve the matter 
overdensity for a given f(T) model. 



III. DYNAMICAL DARK ENERGY AND ITS 
f(T) EQUIVALENT 

In this work we are interested in constructing the f(T) 
scenario that effectively leads to the same background be- 
havior with a given cosmological model. We stress that 
this procedure has a full generality, that is for any dy- 
namical dark energy (DDE) scenario we can construct its 
f(T) equivalent, as long as we know the evolution of the 
dark energy density pdde- However, despite the coinci- 
dence of the background behavior, the perturbations can 
be used to distinguish the given DDE model from the 
corresponding f(T) reconstruction. 



A. Reconstructing the corresponding f(T) family 
for any given Dynamical Dark Energy 

In order to perform the aforementioned reconstruction 
at the background level we start from the pdde effective 
definition flU). Using the fact that / = -YlfHH (since 
T = —6H 2 ), we can extract a differential equation for / 



in terms of Pdde, namely: 
H 



f 



H 



(/ + f 67tGpdde) 



(31) 



The solution to this first-order differential equation can 
be written in closed form as 



f(H) = 1%-kGH J ^j^dH + CH, 



(32) 



where C is an integration constant, and the correspond- 
ing f(T) can be straightforwardly obtained substituting 

describes a 



H by y— T/6. In summary, relation 
one-parameter family of solutions, characterized by the 
parameter C, which by construction mimics perfectly the 
background evolution of the given dynamical dark en- 
ergy density. A self-consistency test for this is that if 
one use the above f(T) and the given H to calculate 
the equation-of-state parameter for the f(T) equivalent 
through (|T5jl . he finds exactly the same result as using 
the conservation equation Pdde = — 3H pdde{1 + w) to 
calculate the equation-of-state parameter for the given 
DDE model. Finally, for the special case of General Rel- 
ativity with a cosmological constant (A) (where pdde = 
PA = A/(8nG)), with C = 0, the corresponding f(T) 
model reduces to the expected / = — 2A. 

Let us now make a comment on the parameter C . In- 
terestingly, and as expected, the C term in (|32p . which 
corresponds to a term proportional to (— T) 1 / 2 in /(T), 
does not have any effect on the background dynamics, 
since it disappears from both Friedmann equations (jlip 
and (fT"2")) . and consequently from any other background- 
level quantity such as luminosity distances or the equa- 
tion of state. In other words, the various specific scenar- 
ios of the family of solutions given by (|32|) . are indistin- 
guishable from each other at the background level. 



B. Using perturbations to distinguish between 
Dynamical Dark Energy and the corresponding f(T) 
family 

In the previous subsection, we showed that the DDE 
and the corresponding f(T) models are perfectly indis- 
tinguishable at the background level. However, this is 
not anymore the case if one takes into account the per- 
turbations, as we show in this subsection. In particular, 
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the perturbations can break the above degeneracy, that 
is the growth history is different for the DDE as well as 
for each member of the corresponding reconstructed /(T) 
family. 

Let us now proceed to the investigation of the pertur- 
bations. For convenience, in the following we express the 
dimensional constant C of relation (|32p in terms of the 
dimensionless quantity 



Cm — — 



CHi 



167rGpDDE,j 



(33) 



where Hi is the Hubble parameter at an initial redshift Zi, 
that is Hi = H(z = Zi), where we use the redshift z as the 
independent variable instead of the scale factor (1 + z = 
cio/ a with ao the present scale- factor value), pdde.i is the 
initial energy density of the DDE. If we assume pdde,« 
to be slowly varying during matter domination, then the 
value of the integral in (|32p at the initial time (z = Zi) can 
be approximated by — 167tGpdde,o an d then Cm can be 
physically interpreted from the fact that the initial value 
of f{H) at z = Zi becomes — 16itG(Cm + 1)pdde,i- 

Amongst the individual members of a reconstructed 
family of f(T) models, corresponding to a given DDE 
scenario, the growth histories can be different for differ- 
ent choices of Cm- In particular, the deviations in the 
growth history for different values of Cm are larger for 
smaller scales, while they diminish for large scales. This 
behavior can be qualitatively understood from equation 
(|30p . Using Sf to denote the matter perturbations under 
the reconstructed f(T) scenarios, we perform in (|3U|) the 
transformation f(T) : f(T) -> F(T) + CH, (where F(T) 
is clearly the Cm = member of the family). We find 
that (|3"D")) transforms to 

3H (l + F' - 12H 2 F") fc 
+ [(3H 2 + k 2 /a 2 ) (1 + F') - 36H 4 F"] cp k 
C k 2 



4:irGp m 6f k 



12H a 



;<t>k- 



(34) 



As we observe, the right-hand-side of the above equa- 
tion indicates the difference between the Sf of the Cm = 
model and that of the various Cm ones. It is there- 
fore clear that for large scales (k —> 0) the difference 
between the Cm and the Cm = models becomes 
small. In summary, perturbations indeed uniquely dis- 
tinguish a DDE scenario from the different members of 
the corresponding reconstructed /(T) family. 

Finally, we mention that our results are similar in spirit 
to [49] in the context of f(R) gravity, where the authors 
find that any expansion history of the universe can be 
replicated by an one-parameter family of f(R) models, 
characterized by a parameter B oc f"(R). which provides 
a variety of different behaviors for the different members 
of the family. However, in the case of the present work 
there is a difference, namely the various models of the 
family have a relatively similar evolution even for large 
differences in the parameter Cm- 



IV. APPLICATION: f(T) EQUIVALENT FOR 
QUINTESSENCE 

In the previous section we described how we can con- 
struct f{T) models that exhibit the same behavior at 
the background level with any given dynamical dark en- 
ergy (DDE) scenario. Furthermore, we showed how the 
examination of perturbations can be used in order to dis- 
tinguish the DDE scenario from the various members of 
the reconstructed f(T) family, as well as these various 
members from each other. In the present section, in or- 
der to provide specific examples, we apply these proce- 
dures in the well-known dynamical dark energy scenario 
of quintessence. 

In the quintessence paradigm, the dark energy sector 
is attributed to a homogeneous scalar field $, while the 
matter sector is described by a pressureless perfect fluid 
with energy density p m [13l - fl6| . In the case of flat geom- 
etry the Friedman equation reads 



H = 



8ttG 



$2 

Pm + -S- + V(*) 



(35) 



where U($) is the scalar-field potential. In these mod- 
els, the field typically rolls down a very shallow potential, 
eventually coming to rest when it can find a local mini- 
mum. In particular, the evolution equation for $ reads: 



$ + 3i/$ + U'($) = 0, 
while the corresponding one for p m is 
Pm + 3Hp m = 0. 



(36) 



(37) 



Additionally, one can define the effective dark energy 
density and pressure in field terms as 



Pdde 



Pdde 



$2 

~2~ 



(38) 
(39) 



and thus the corresponding dark-energy equation-of-state 
parameter as 



WDDE 



Pdde 
Pdde 



^ + u($) : 



(40) 



which obviously present a dynamical behavior in general. 
Finally, in such a scenario, the growth of perturbations 
is g overned by the linearized Einstein equations (see e.g. 



v f 




5$ = -3H6& - 



V") 6$ + 4$0 fc - 2V'4>k, (41) 
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where Vf is the velocity potential of the matter fluid and 
primes denote derivatives with respect to 4>. 

One last subject that has to be settled before proceed- 
ing forward is the choice of the quintessence potential 
V($>). Amongst the various ansatzes of the literature, 
in this work we focus on three commonly studied cases, 
namely: 

1. The Pseudo-Nambu-Goldstone-Boson (PNGB) 
Model: 

V ($) = A 4 [1 + cos ($//)] , (42) 

where A is the energy scale of the potential, and / is 
a symmetry restoration scale. This model was first 
proposed in (5l| , while its cosmological applications 
were discussed in (H3 |. 

2. The Exponential Model: 

V ($) = A 4 cxp (-0<S>/M P1 ) , (43) 

where Mpj is the Planck mass. Exponential poten- 
tials arise in a variety of contexts, such as higher di- 
mensional gravitational theories, moduli fields, and 
also non-perturbative effects such as gaugino con- 
densation (see e.g. [53j and references therein). 

3. The Power Law Model: 

V (*) = A 4+n <S>- n . (44) 

Power law potentials have been studied extensively 
in fljjl . 0, 0] , and have been shown to arise in the 
context of SUSY in 

The usual method of differentiating various types of 
DDE models is to examine their phase-space behavior, fo- 
cusing in particular on the dark-energy equation-of-state 
parameter and its evolution (56| . The various models can 
either fall into a "thawing" or "freezing" category. The 
former is thusly named because the field is only recently 
beginning to exhibit dynamical behavior since it is being 
frozen at some distance from its minimum due to Hubble 
friction, while the name of the latter arises form the fact 
that it is dynamical throughout much of the universe's 
history but it becomes frozen into place at late times 
(during the time of dark energy domination). Amongst 
the three choices above, the first two are thawing models 
while the third one is a freezing model. 

A. Reconstructing f(T) for quintessence 

Let us now reconstruct an f(T) cosmology that ex- 
hibits the same behavior at the background level with a 
given quintessence scenario. As described in subsection 
IIH Al the corresponding f(T) family of models is given 
by relation (|3"2"|) , where in the quintessence case the dark 
energy density pdde is given by (|38p . while the cosmo- 
logical equations close by (|Bl))) - (f3"T|) . 



In a non-trivial quintessence scenario, (|3"2"T) cannot 
be solved analytically. Thus, in the following we per- 
form a numerical elaboration in the case of the three 
quintessence models described above. In particular, 
knowing the evolution of pdde and of the effective 
dark-energy equation-of-state parameter w for the three 
quintessence models, we numerically reconstruct the cor- 
responding f(T) family using (|32[) . and then we numeri- 
cally extract the w-behavior using (|T51) . 

For the three quintessence potentials we suitably 
choose the parameter values in order for w to lie within 
its observational limits -1<uidde<-0.9 for low redshifts 
[H . We solve the cosmological equations determining the 
initial conditions at a redshift Zi deep inside the mat- 
ter dominated era, and we evolve the system until the 
density parameter of matter il m (z) = 8irGp m /(3H 2 (z)) 
becomes equal to ss 0.3 at z — 0, as it is required by ob- 
servations. For our numerical analysis, we work in units 
where 8irG = 1 and p\ = 4/3. 
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Figure 1: For members of the reconstructed f(T) family 
of models for the Pseudo-Nambu-Goldstone-Boson (PNGB) 
quintessence scenario given by \4ty , characterized by four 
choices of the parameter Cm . 

In Figures [T][3] we present the reconstructed f(T) 
for the three quintessence models, namely the Pseudo- 
Nambu-Goldstone-Boson (|4"2"T) . the exponential (|43p . and 
the power-law (|44[) models, respectively (such figures are 
easily created since we know numerically f(z) and H(z), 
that is T(z), and thus we obtain f(T)). In each figure, 
we depict f(T) for four members of the f(T) family of 
models, characterized by four choices of the parameter 
Cm- 

We mention that the C-term in (|32| behaves like a de- 
caying mode, as it falls of approximately as a~ 3 / 2 during 
early times, making the C = (Cm = 0) case an attrac- 
tor for this family of models. This is demonstrated in 
Figures [T][3l where it is clear that the Cm models 
approach the Cm = one, for the reconstructions of all 
three quintessence scenarios. Finally, note that we have 
focused only on positive CM-values, since negative val- 
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Figure 2: For members of the reconstructed f(T) family 
of models for the exponential quintessence scenario given by 
\4ty , characterized by four choices of the parameter Cm ■ 
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Figure 4: The evolution of the dark-energy equation-of- 
state parameter w(z), as a function of the redshift, for the 
Pseudo-Nambu-Goldstone-Boson (PNGB) quintessence sce- 
nario given by as well as for four members of its corre- 
sponding reconstructed f(T) family of models (characterized 
by four choices of the parameter Cm ). 




Figure 3: For members of the reconstructed f(T) family of 
models for the power-law quintessence scenario given by \Jjl$ , 
characterized by four choices of the parameter Cm . 



ues will produce graphs which are symmetrically reflected 
about the f(z) = f (z = 0) = -2A line. 

Having reconstructed the f(T) family of models that 
corresponds to the three quintessence scenarios, in Fig- 
ures HUH we present the behavior of w(z), which is a ba- 
sic observable, for each potential respectively, according 
to (|4l?)) . In each figure, we depict additionally the w(z) 
of the corresponding reconstructed f(T) family of mod- 
els, according to (fT5|) . for four values of the parameter 
Cm- As we observe, in all cases there is a perfect overlap 
of the corresponding figures, which verifies the fact that 
for background-level quantities, such is the dark-energy 
equation-of-state parameter w, any dynamical dark en- 
ergy scenario is indistinguishable by construction from 
its corresponding f(T) equivalent family of models. 
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Figure 5: The evolution of the dark-energy equation-of-state 
parameter w(z), as a function of the redshift, for the expo- 
nential quintessence scenario given by [43ty , as well as for 
four members of its corresponding reconstructed f(T) fam- 
ily of models ( characterized by four choices of the parameter 
Cm )■ 



B. Using perturbations to distinguish between 
quintessence and the corresponding f(T) family 

In the previous subsection, we applied our general 
/(T)-reconstruction formalism of section lTlI Al in the case 
of quintessence dynamical dark energy scenario. In the 
present subsection we investigate the perturbations fol- 
lowing subsection IIIIBl in order to distinguish between 
quintessence and its equivalent family of models. 

The perturbation evolution for quintessence is given 
by equations (|4"Tj) . while for the family of f{T) models 
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Figure 6: The evolution of the dark-energy equation- of- state 
parameter w(z), as a function of the redshift, for the power- 
law quintessence scenario given by \44\ )> as well as for four 
members of its corresponding reconstructed f(T) family of 
models (characterized by four choices of the parameter Cm )■ 



by (|34p . We follow the numerical elaboration described 
in the background case, and we additionally choose the 
initial velocity of the metric perturbation to be at zero. 



Cm = 0). Similarly, in Figures \§\ and [TUl we provide the 
corresponding plots for the exponential scenario given by 
(|4U)) . while in Figures [Til and we show the correspond- 
ing graphs for the power- law scenario given by (1441) . 
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Figure 8: The evolution of the matter overdensity S as a 
function of the redshift z, on a scale of k = lh Mpc -1 , for 
the PNGB scenario given by as well as for three mem- 

bers of its corresponding reconstructed f(T) family of models 
(characterized by three choices of the parameter Cm )■ 
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Figure 7: The evolution of the matter overdensity 8 as a 
function of the redshift z, on a scale of k — 10 _3 /i Mpc -1 , for 
the PNGB scenario given by \4ty , as we ^ as f or three mem- 
bers of its corresponding reconstructed f(T) family of models 
(characterized by three choices of the parameter Cm)- 

In Fig. [7] we demonstrate the evolution of the matter 
overdensity 5 as a function of the redshift, in the case of 
Pseudo-Nambu-Goldstone-Boson (PNGB) quintessence 
scenario given by (|42|) . in a scale of k = 10~ 3 h Mpc -1 , 
while in Fig. [5] we provide the same graph but in a scale 
of k = lh Mpc . In the same Figures, we additionally 
depict the (5-evolution for three members of the corre- 
sponding reconstructed /(T) family, determined by three 
distinct values of the parameter Cm {Cm — ±5000 and 
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Figure 9: The evolution of the matter overdensity 5 as a func- 
tion of the redshift z, on a scale of k = 10~ 3 h Mpc -1 , for the 
exponential scenario given by \4$il , is well as for four mem- 
bers of its corresponding reconstructed f(T) family of models 
(characterized by four choices of the parameter Cm )■ 



From these Figures we observe that the evolution at 
the perturbation level distinguishes the quintessence sce- 
nario from its reconstructed equivalent family of /(T) 
models, as well as the various f(T) family members from 
each other. Concerning the scale, as discussed in subsec- 
tion IIII Bl the deviation in 5-evolution is strong for small 
scales and weak for large scales. Concerning the redshift, 
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Figure 10: The evolution of the matter overdensity 5 as a 
function of the redshift z, on a scale of k = h Mpc^ 1 , for the 
exponential scenario given by \4'ty , as well as for four mem- 
bers of its corresponding reconstructed f(T) family of models 
(characterized by four choices of the parameter Cm ). 



Figure 12: The evolution of the matter overdensity S as a 
function of the redshift z, on a scale of k — h Mpc^ 1 , for the 
power-law scenario given by \44ty i as well as for three mem- 
bers of its corresponding reconstructed f(T) family of models 
(characterized by three choices of the parameter Cm )■ 



to 7 





— 5 




<--s f rc M =o] 




— S f [C M =5000] 




*P*F sooo; 







1 1.5 2 2.5 3 

1+z 



Figure 11: The evolution of the matter overdensity 8 as a 
function of the redshift z, on a scale of k = 10~ 3 h Mpc^ 1 , 
for the power-law scenario given by \44% as well as for three 
members of its corresponding reconstructed f(T) family of 
models (characterized by three choices of the parameter Cm)- 



for all the models the deviation is negligible for large red- 
shifts z > 3, which is expected since the universe is mat- 
ter dominated at that time and the dark-energy sector 
(DDE or f(T)) plays no role in determining the growth 
of perturbations. Moreover, note that that the shape 
of the (5(z)-curve differs fundamentally from that of the 
5f(z; CM)-curves. Therefore, it is possible to choose a 
value of Cm such that S and Sf coincide at a particular 
redshift, but they will not coincide at other redshifts for 
this particular choice of Cm- 

Finally, as we can see, in general the parameter Cm has 
a very weak impact on the splitting of 5-evolutions, even 
for small scales - the difference between the Cm = 5000 



and the Cm = —5000 cases is modest in all the cases 
considered. This is most likely a result of the attractor 
nature of the Cm = solution described in IIV A[ 



V. CONCLUSIONS 

In this work we investigated f(T) cosmology in both 
the background, as well as in the perturbation level, as 
a way to mimic the behavior of a given model of dy- 
namical dark energy. f(T) gravity is the extension of the 
"teleparallel" equivalent of General Relativity, which uses 
the zero curvature Weitzenbock connection instead of the 
torsionless Levi-Civita connection, in the same lines as 
f(R) gravity is the extension of standard General Rela- 
tivity. 

First, we presented the general formalism for recon- 
structing the f(T) scenario that effectively leads to the 
same behavior with any given cosmological model, at the 
background level. As we showed, for a given dynamical 
dark energy scenario, one can reconstruct its equivalent, 
one-parameter family of models. Although they exhibit 
completely indistinguishable background behavior, the 
perturbations can break the above degeneracy. That is, 
the growth history is different for the given cosmological 
scenario as well as for each member of the corresponding 
reconstructed f(T) family. 

In order to present our results more transparently, we 
applied the aforementioned general formalism in the well 
known case of quintessence dynamical dark energy sce- 
nario. After reconstructing the /(T) family for three 
quintessence scenarios, we numerically extracted the evo- 
lution of the dark-energy equation-of-state parameter, 
which is a basic observable. Since it is a background 
quantity it leads to identical behavior for quintessence as 
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well as for its f(T) equivalent. 

However, upon examining the perturbations, and in 
particular the evolution of the matter overdensity, we do 
acquire a way to distinguish between quintessence and 
its f(T) equivalent family, as well as between the infi- 
nite members of the reconstructed f(T) family. More 
specifically, the deviation is strong for small scales and 
weak for large scales, and additionally it is negligible for 
large redshifts (z > 3), since there the universe is matter 



dominated with a negligible dark-energy sector. 

In summary, f(T) gravity can effectively mimic any 
dynamical dark energy cosmological scenario at the back- 
ground level. The perturbation analysis can break this 
degeneracy leading to rejection or acceptance of specific 
f(T) models. These features make f(T) cosmology an 
interesting candidate for the description of nature, which 
requires deeper examination. 
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